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In an organic blend the vibrational normal mode excited by exciton splitting is the same as the 
one coupled to charge hopping. Excess driving force for exciton splitting can therefore aid charge 
transfer, if vibrational relaxation is slow compared to charge transfer. A model is developed that 
takes this into account and hence explains the experimentally observed relation of driving force for 
exciton splitting and charge yield and that high charge yields can be achieved with the observed 
fast rates of recombination. 
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INTRODUCTION 

In an organic blend, charges are generated by split- 
ting excitons at the interface between materials with dif- 
ferent electron affinities: the donor material is typically 
responsible for absorbing light and has a greater affin- 
ity for holes; the acceptor has greater affinity for elec- 
trons. Charges must escape each other's Coulomb at- 
traction without recombining. Exciton splitting is typi- 
cally driven by a large (> 0.5 eV) driving force We 
develop a model where this excess energy aids charges 
escape from the Coulomb well. This can happen because 
the polaron pair formed by exciton splitting is in a vi- 
brational excited state. The mode excited is the super- 
position of the modes that change the donor from the 
exciton to the charged geometry and the acceptor from 
the neutral to the charged geometry. This second mode is 
also coupled to charge transfer in the acceptor molecule, 
hence, if the mode is sufficiently long lived, the states 
formed just after exciton splitting are able to separate 
more efficiently if extra energy is available. 

That hot states are involved in charge separation is not 
a new idea. Peumans et al. need to assume very large 
thermalization radii to reproduce experimental charge 
yields Q. Models involving hot excitons have been used 
to describe charge generation in pure polymers Q . Mod- 
els by Offermans et al. predict increased charge yields 
if charges are generated hot Experiments have also 
shown that charge mobilities are higher immediately af- 
ter charge generation Q , suggesting that charge transfer 
is fastest immediately after exciton splitting. However, 
how thermalization lengths are connected to molecular 
properties is not yet understood. 

Understanding the link between the driving force for 
charge separation and quantum yield will help determine 
the limiting efficiency of organic solar cells (OSC). OSCs 
can have high quantum efficiencies Q , however open cir- 
cuit voltages (Voc) are still substantially smaller than the 
band gap. The Voc for several materials is less than half 
the maximum of external quantum efficiencies Until 
recently, designing better OSC materials was thought to 
require reducing the energy lost in exciton splitting Q. 



Recent measurements by Okhita et al. [8[ relate the en- 
ergy loss and the yield of charges: bad news for OSC as 
increasing the Voc would also lower quantum yields. 

Modeling generation without thermalization requires 
recombination rates in the microsecond range m mm 
in order to achieve high yields of charges (> 10%). But 
experimental measurements of the rate of charge recom- 
bination have found it to be faster than nanoseconds 
12 1 13 • Here high yields are obtained with fast recom- 
bination by including the effect of charge transfer from 
hot states. 

The theoretical background for the model developed 
is given by the Marcus Levitch Jortner (MLJ) equation 
[14j . The rate of charge transfer becomes: 
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where (i\H\j) is the matrix element for charge tunnel- 
ing between the diabatic electronic states i and j, X is 
the outer sphere reorganization energy four coupling to 
a classical bath, uj is the frequency of the quantum mode 
and (0|/) is the Franck Condon overlap between the ini- 
tial vibrational state (the ground state) and the final one 
(with quantum number f). The Franck Condon factor is 
computed assuming that the initial and final vibrational 
states are wavefunction of simple harmonic oscillators 
with identical frequencies but shifted in mass weighted 
coordinates by y/A, the Huang Rhys factor. AG rep- 
resents the difference in free energy between the initial 
and final diabatic states. Each term in the summation 
represents the probability of resonant electron tunneling 
occurring between the ground vibrational state of the ini- 
tial state and the / vibrational state of the final states. 
This equation can be used to model all three processes 
of interest in the paper: exciton splitting, charge trans- 
fer and charge recombination. A caveat is that charge 
transfer must occur in the non-adiabatic regime, i.e. the 
transfer integral must be small compared to the electron 
A Q and to Atku. 

Figure [T] shows the processes modeled in the paper. 
Each column of energy levels represents an electronic 
state, labeled by the distance d between electron and 
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hole. At a distance the system is in the neutral ground 
state (black energy levels) or in the excited state (red 
energy level). The processes competing are: i) exciton 
dissociation; ii) charge transfer from a vibrationally ex- 
cited state; iii) decay of the vibrational energy by internal 
conversion; and iv) charge recombination to the ground 
state. Exciton splitting is not modeled explicitly, but is 
used to provide the initial excited vibrational state. 



METHOD 

The vibronic state \i, v) of the system is described by 
the outer product of the electronic state \i) and of the 
vector describing the vibrational state of the oscillators 
\v). Each mode in \v) is coupled to the energy of a single 
electronic state 

When charge transfer occurs between two electronic 
states i and j , only the modes of the oscillators coupled 
to each of those two states will be able to change. If 
the mode associated with electronic state % changes from 
quantum number p to p and the mode associated to j 
changes from v to a, the rate of charge transfer can be 
written as: 

T ( i,p,,v -> j,p,cr) = (2) 
]T P(S, A\v, p)T(i, S -> j, S')P(S', A\p, a) 

S'.S 

where T(i, S —> j,S') represents the rate of transfer be- 
tween electronic states i and j and symmetric combina- 
tion of the vibrational modes S to S'. Only symmet- 
ric combinations of modes contribute to charge transfer, 
as the energy difference between two oscillators depends 
only on the difference of the distorsions of each oscillator. 
The probability terms P{S,A\p,v) represent the proba- 
bility of a system of two simple harmonic oscillators in 
localised modes p and v to be in a symmetric mode S and 
an antisymmetric mode A. By using a probability term, 
rather than explicitly considering the wavefunction, we 
implicitly assume that with each charge hop, coherence 
is lost. 

The rate term T(i, S — > j, S') can be computed with a 
trivial generalisation of the MLJ equation: 
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where all symbols have the same meaning as in equation 
1. 

The probabilities P(S, A\p, v) can be determined by 
assuming that the ladder operators for the symmet- 
ric/antisymmetric modes are given by combinations of 
the appropriate localised ladder operators a\ and a\: 
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where the 5 functions ensures that the probability is non- 
zero only if p + v = A + S. 

As well as charge hopping between vibronic states, a 
state is also allowed to lose a vibrational quantum by 
internal conversion at a rate kyR without changing its 
electronic state. 

The dynamical system defined by these equations is 
linear and could be solved by writing a master equation, 
but since the total vibronic phase space is large it is more 
efficient to solve it by generating a starting state vector 
\\,v) and updating it using a continuous time random 
walk algorithm [15]. Rates are then computed only for 
states as they are needed. The simulation is stopped if 
either the charges recombines, or if they reaches a certain 
distance d. We do not explicitly model exciton dissocia- 
tion, but use the parameters for exciton dissociation to 
determine the probability of starting the simulation from 
a particular vibrational level, given by: 
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where each term in the summation represents the prob- 
ability of producing a certain symmetric state S and the 
probability that the state gives rise to localised state v. 
The excited neutral state is labeled 0*. Note that this 
expression depends on the difference in energy between 
the excited state and the first charged state (labelled AG 
in figure [l}, on A , on the Huang Rhys factor and on fiui, 
but not on the matrix element. 

We apply this model it to a one dimensional chain of 
acceptors, with the hole localised at one end of the chain. 
Electronic states are uniquely labeled by the distance be- 
tween hole and electron \i), for example state |1) will have 
the electron and hole one lattice apart spacing. State |0) 
is the neutral ground state reached if the electron and 
hole recombine. Each state \i) is connected to its neigh- 
bors \i + 1) and \i — 1) with the same transfer integral 
for charge separation Vcs, with the exception of states 
|1) and |0), which are connected by the transfer integral 
for charge recombination Vcr, this allows the timescales 
for charge transfer and recombination to be controlled 
independently. 

The energy E(i) of a particular electronic state \i) is 
determined only by the applied electric field F and the 
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FIG. 1: Sketch of the mechanisms described in the main text: 
i) exciton splitting, ii) charge transfer from a vibrational ex- 
cited state, iii) internal vibrational relaxation, iv) recombi- 
nation. The red arrow is a trajectory leading to charge gen- 
eration whereas the blue one to recombination. The pseudo 
bandgap Et and the driving energy for exciton splitting AG 
are shown. 



Coulomb potential: E(i) 



F a i, where 
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a is the lattice constant, e and eo are respectively the 
relative and vacuum permittivities and e is the charge 
of an electron. The energy of state |0) is a fixed value 
Eb below the zero of this potential. Ej, is the the pseudo 
band gap between the highest occupied molecular orbital 
of the donor and the lowest unoccupied molecular orbital 
of the acceptor. Eb and AG (the driving force for exciton 
splitting) are shown schematically in figure [TJ 



RESULTS 

All our modeling uses the following parameters for the 
energetics: a = 1 nm, Eb = —0.9 eV, e — 4 and F = 
5 10 5 V/ cm. The electric field employed is rather large, 
smaller fields would be necessary if more subtle models 

1(5. 17] lowered the 



of the polarization of the interface 
barrier do charge separation, or if a three dimensional 
model which allowed both charges to move was employed 

The parameters of the charge transfer equation used 
are: A = 1, hu = 0.17 eV, A c = 0.1 eV, V CS = 0.05 eV, 
and kvR — 10 12 s -1 , a set of values typical for fullerene 
[l8j ]. The parameters for recombination are the same as 
above, but with a transfer integral Vcr = 0.01 eV and 
an energy defined by Eb = — 1 eV. This puts the recom- 
bination rate in the sub ps regime. Panel a) in figure [5] 
shows the conditional probability of the system reaching 
a certain distance d given that it has already reached a 
distance d — 1 for a range of initial vibrational quanta. 



Showing the conditional probability helps identify the 
thermalization distance, because once the system ther- 
malizes all the curves fall onto each other. Increasing the 
number of phonons in the initial state massively increases 
the probability that the charge is able to escape recom- 
bination, for example in the vibrational ground state the 
probability of rccombining immediately is 60%, but with 
just one quantum of vibrational energy this probability 
decreases to less than 1%. Simulations with one extra 
phonon thermalise after just one hop. The higher excited 
states thermalise at a distance of 4nm, showing that (for 
these parameters) thermalization is complete after three 
or four hops. The inset of panel a) shows the probabil- 
ity of a particular vibrational excited state being created 
given a certain AG . Clearly highly excited states are 
more likely for greater AG. Weighing the probabilities 
computed by the probability that that starting number 
of phonons is achieved for a certain AG gives the charge 
yield as a function of AG. 

Panels b),c) and d) show the effect on charge yield 
at 5 nm as a function of AG of: the transfer integrals 
for charge transport (b), the rate of vibrational relax- 
ation (c) and A D for charge transport (d). Loss of vi- 
brational excitation can occur in two ways: by internal 
vibrational relaxation or by losing energy to the clas- 
sical modes through the external reorganization energy. 
The yield of charges is therefore the result of competi- 
tion of three basic process: charge transport, internal 
vibrational relaxation and charge recombination. De- 
creasing the transfer integral not only reduces the yield of 
charges, it also decreases the dependence of yield on AG. 
This is because slowing down charge transfer not only 
makes recombination more likely (reducing the yield) it 
also makes vibrational relaxation faster (reducing the de- 
pendence on AG). Increasing the rate of vibrational re- 
laxation reduces the yield because it causes relaxation 
to occur faster. A has a large effect because it dictates 
how much energy is lost to the classical bath with each 
hop: making it larger makes the system thermalise more 
effectively through charge transfer. 



CONCLUSION 

Explicitly modeling charge transfer between vibra- 
tionally excited states allows high charge yields with fast 
recombination and explains the dependence of charge 
yield on the driving force for exciton splitting. Improv- 
ing charge transport improves charge generation because 
of more effective competition of charge separation with 
both charge recombination and vibrational relaxation. 

Three conditions must be fulfilled for vibrational exci- 
tation to aid charge generation : 1) the vibrational mode 
excited by exciton splitting must be the same as that 
coupled to charge separation, 2) vibrational relaxation 
must be slow compared to charge tunneling, 3) A D must 
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FIG. 2: Panel a) The conditional probability of the system reaching a certain distance d given that it has reached distance 
d — 1 for different initial quanta of vibrational excitation. The inset shows the probability distribution of producing a certain 
number of quanta for various values of the driving force for exciton splitting. The curves are offset for clarity. The value of 
AG in eV is shown below each curve. Panels b), c) and d) show the yield of charges at 5 nm as a function of the following 
parameters: b) the transfer integral for charge separation, c) the rate of internal relaxation, d) the outer sphere reorganization 
energy. Each curve is labelled by the value of the parameter used in that simulation. 



be small to ensure low losses to the classical bath. 

The present model makes two fundamental by two as- 
sumptions: that coherence is lost at each hop and that 
the coupling is non-adiabatic. Future work concentrate 
on tackling these problems. The model allows us to pre- 
dict that the quantum efficiency of a blend should in- 
crease with photon energy. We have assumed that vibra- 
tional excitation aids charge separation only in the accep- 
tor molecule. If the modes that lead to exciton splitting 
are not orthogonal to those coupled to charge transfer 
also in the donor molecule, the effect of the initial energy 
would be increased. This provides a design rule for good 
donor materials. 
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